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Abstract 

The eight-port homodyne detection apparatus is analyzed in the framework 
of the operational theory of quantum measurement. For an arbitrary quan- 
tum noise leaking through the unused port of the beam splitter, the positive 
operator valued measure and the corresponding operational homodyne ob- 
servables are derived. It is shown that such an eight-port homodyne device 
can be used to construct the operational quantum trigonometry of an optical 
field. The quantum trigonometry and the corresponding phase space Wigner 
functions are derived for a signal field probed by a classical local oscillator 
and a squeezed vacuum in the unused port. 



PACS number(s): 42.50.Dv, 42.50.Ar 



1 



Typeset using REVTpjX 



I. INTRODUCTION 



In this paper we present a general operational description of measurements performed 
with the help of the eight-port homodyne detector. This experimental setup has been 
extensively used in the past years in connection with the homodyne detection of optical 
signals |l| and the determination of the phase of quantum states of light [0. The series 
of beautiful experiments, performed by Noh, Fougeres and Mandel (NFM) has provided a 
fully operational approach to the measurement of the quantum phase of a single mode light 
field. These experiments have inspired a lot of further theoretical works devoted to a general 
operational approach to quantum measurements in quantum optics. 

The NFM experiments have been treated theoretically in many different ways 
however our approach will be entirely based on the concept of the operational phase space 
propensity [[F] and on the associated with it operational quantum observables ||. Such a 
formulation, which is one of many possible descriptions of realistic measurements, is closely 
connected with the concept of the positive operator valued measure (POVM). In the op- 
erational formulation the POVM replaces the spectral measure of the intrinsic quantum 
observable PJlO|]. 

The operational formalism has been applied to the NFM experimental device in the 
context of the operational definition of the Hermitian phase operator |J. In the framework 
of such an approach, the so called quantum trigonometries for the optical fields have been 
constructed [|TT|-|T3|]. Recently the operational theory has been applied to the homodyne 
detection scheme with a fixed or a random phase of the local classical oscillator |14]] . 

The advantage of using the operational measurement theory is that we can easily and 
naturally take into account the influence of the experimental device on the measured system. 
For example, the eight port homodyne detector is characterized by the efficiency of the 
detectors, the reflectance and the transmittance of the beam splitters, and the state of the 
quantum field entering the unused ports of the device. This additional "noise" field plays 
a fundamental role, if the NFM device is used to investigate the phase of quantum optical 
fields. In such a setup the impact of the correlations on the probe field and the noise field 
plays a crucial role in the operational definition of the quantum phase. This particular 
interplay between the quantum signal and the quantum noise, leaking through the unused 
port in the NFM device, provides the framework of the operational approach to the quantum 
phase. The case of the noise field being the vacuum state has been analyzed, in the spirit 



of the operational formalism, in Ref. |11J| . These results have been generalized for the case 
of the squeezed vacuum state In this paper we give a more detailed discussion of the 
squeezed vacuum state case and present general results obtained for an arbitrary noise field. 

In Section II we shortly review the operational approach to the measurement in quantum 
mechanics. In Section III we derive the operational probability density — propensity for 
an arbitrary noise field in the NFM apparatus. In Section IV we derive the operational 
homodyne observables for the eight-port homodyne detection scheme. Section V is devoted 
to the discussion of the quantum trigonometry, generated in the operational phase measure- 
ments in the NFM device, with the noise being a squeezed vacuum. Finally, some concluding 
remarks are presented. 
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II. OPERATIONAL APPROACH TO MEASUREMENTS IN QUANTUM 

MECHANICS 



The approach to the quantum mechanical theory of measurement that will be used in 
this paper was described in details in Ref. |§. Below we just recall the main results of 
this work in order to make the paper self-contained. The statistical outcomes of an ideal 
measurement of a certain observable A: A\a) = a\a) are described by the spectral measure 

P 4a) = \(am 2 , (1) 

where l^) G 7i is the state vector of the measured system. It is known that the spectral 
measure contains all the relevant statistical information about the investigated system but 
it makes no reference to the apparatus employed in the actual measurement. Due to this 
property A will be called an intrinsic quantum observable. 

The description of more realistic measurements necessarily involves additional degrees 
of freedom such as, for example, the unused ports of the beam splitters or the local os- 
cillator being the phase calibrating device. In such an approach, in order to describe the 
quantum mechanical measurement, we have to work in the extended Hilbert space 7i ® 'Hp, 
containing both the Hilbert space of the investigated system 7i and the Hilbert space Tip of 
all the additional measuring devices, which will be called for short a filter. The process of 
measurement is described by the interaction between the system and the filter. Because we 
are interested only in the properties of the system, we shall reduce the degrees of freedom 
by tracing the dynamical evolution in the combined Hilbert space over the filter Hilbert 
space obtaining in such a way an operational propensity Pr(a), which is nothing else but a 
probability distribution of a certain classical quantity a measured in a real experiment. This 
propensity, for arbitrary density operator p of the system, is 

Pr(a) = Tr{^(a)p}, (2) 

with the normalization condition 

J da Pr(a) = 1 J dafi(a) = 1, (3) 

where fi(a) is a filter dependent POVM. We see that in a realistic measurement, involving 
a filter, the spectral decomposition of A is effectively replaced by a positive valued measure 
daf{a) §. 

In view of the linear relation between the propensity and the POVM, the operational 
statistical moments of the measured quantity 

= J daa n Pr(a), (4) 

define uniquely an algebra of operational operators 

^=<4 n) >, (5) 

where 
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AP = J daa n f(a). (6) 

These operators are called operational operators, because they represent quantities measured 
in a real experiment i.e., described by a dynamical coupling of the measured system and the 
filter. As a rule, the algebraic properties of the set A^ differ significantly from those of the 
powers of A. For instance, (A^p) 2 ^ . This property will have important consequences 
in the discussion of the uncertainty introduced by the measurement. 

It is seen that the propensity and the operational operators are natural equivalents of 
the spectral probability distribution and the intrinsic operators. The difference being that 
the operational observables carry information about the system under investigation and the 
selected measuring device, which is represented in this paper by the eight-port homodyne 
detector. 



III. NFM APPARATUS WITH AN ARBITRARY NOISE FIELD 

In this Section we derive a general formula for the POVM that describe the NFM device 
with an arbitrary noise filed leaking through the unused port of the beam splitter. We shall 
assume that the NFM apparatus has unit detectors efficiency, that a 50%: 50% beam splitters 
are used and that the local oscillator is classical. We assume that we want to measure an 
observable A(b, w). In general such an observable requires additional informations about the 
particular ordering of the boson creation and annihilation b operators. For example the 
quantum observable that corresponds to the polar decomposition of the amplitude of a field 

mode could be defined as bj \fb^i), however this formula is not unique because of the possible 
different ordering of the boson operators. Due to the ordering problem, there is no canonical 
way of selecting the right, i.e., the physically justified ordering. The operational approach 
removes this ordering ambiguity. We shall perform a measurement of this observable using 
the NMF device. In this case the filter Hilbert space is reduced to a Hilbert space 7i v , 
where v denotes the quantized mode leaking through the unused port. In the extended 
filter-system space TCb <8> TC V the observable A will be the following function of the creation 
and annihilation operators 

A = A{h + v\tf + v), (7) 

where the operators b, represent the probe field and the operators v, v> describe the noise 
field, i.e. the filter degrees of freedom. In this extended space, the combination b + v< do 
commute with W + v and such operators can be measured simultaneously. In the case of the 
unimodular phase operator the combination 

h + " (8) 



is operationally unique and the extended field modes can be simultaneously measurable. A 
similar approach to the measurement of a simple quantum system, using an extension by an 
additional degree of freedom, has been used by Einstein Podolsky and Rosen ||16|| . In their 
case the particle system with the non commuting position and momentum operators p and 
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q has been extended to commuting observables p + P and q — Q, where P and Q can be seen 
as the momentum and the position of the filter. This physical procedure has a counterpart 
in mathematics, and is called the Naimark extension of the POVM into a projective measure 
on a larger Hilbert space |TD . 



The set of operational operators corresponding to the powers of the operator A in the 
extended filter-system are defined as 



AP = Tr v {A n (b + v\ + v)g{v\ v)} . 



(9) 



We will assume that the noise field density matrix g(v', v) can be expressed in the normally 
order form i.e., all the annihilation operators are to the right of the creation operators. 
This is true for all physically interesting cases discussed in the literature. Provided that 
this condition is fulfilled we can rewrite the above expression, using the coherent state 
decomposition of unity, in the following way 



A 



d 2 a 



71 



A n (a, a*)g(a* -b\a-b). 



(10) 



This formula gives us the form of the POVM for the experimental NFM setup with arbitrary 
noise: 



F{ol) = g{a* — ¥,a — b), 



with the normalization: 



d 2 a 



71 



(11) 



(12) 



Such a result has a simple meaning in terms of shifts in the coherent state phase space of 
the signal field. In order to measure the signal mode we have to probe it with a filter state 
shifted in the phase space by the amount a, and sum up over all possible manifolds of shifts 

0- 

As an example of this general approach we find explicitly the POVM for several inter- 
esting cases of the filter state. For the vacuum state 



g(v, v 



|0)(0| 



exp (— v'v 



the corresponding POVM is 



jF = | a) (a | 



where we employed the following identity [ 1I7J 



|a)(at| 



cxp 



-(6 f -a*)(6-a)] : 



(13) 



(14) 



(15) 



and have denoted the normal ordering by semicolons. This result has been obtained already 



in Ref. pT| . Another interesting state is the squeezed vacuum state []rj 

cosh -1 s 



t 



\s)(s\ 



x : exp 



tanh s (t 



e^ t2 + e 



(16) 
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which leads to the POVM in the form 



T = \a, s) (a, s\ 



(17) 



which has been derived in Ref. | I2| . The squeezed state is generated from the ground state 
by the action of the squeezing operator S(s, 0) and the Glauber displacement operator D(a). 
Similarly, for the Fock state |n)(n| we obtain 



T = D(a) 



n + 1 



n 



n 



1 



r I a 



a) 



(18) 



which equals to the POVM from the Ref. 
position and momentum measurement. 



introduced in the context of the operational 



IV. OPERATIONAL HOMODYNE OBSERVABLES AND THE POSITION AND 

MOMENTUM MEASUREMENT 



As a first application of the above described formalism we shall discuss the operational 
measurement of the homodyne quadratures. Such a discussion has already been carried in 
14| , for a homodyne detection in the standard configuration i.e., with one beam splitter 



and two detectors. Here we present the operational observables associated with the joint 
measurements of two field quadratures with the help of the NFM apparatus. We will also 
slightly generalize our previous considerations, assuming that the beam splitter, which mixes 
the probe field with the noise field has an arbitrary transmittance coefficient T and that the 
noise filed is a squeezed vacuum. As before, we shall assume that the local oscillator is a 
strong coherent filed described by a = \a\ exp (i9). 

In the extended filter-system space, the operators leading to the operational observables 
of the two homodyne quadratures are the "unnormalised" sine and cosine operators |||| 



X x = - 
2 



a 



+ a 



and 



Xo 



ia(Vl - TP + VTv) - ia*(y/l-Tb + VTtf)] 



If we define two generating operators (i = 1, 2) for the corresponding quadratures 

Zir(£) = Ti v jexp (£X;) , 
then the operational moments are obtained by the n-fold derivatives: 



X 



(n) 



(19) 



(20) 



(21) 



(22) 



For both quadratures the generating operators can be easily found 
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ZiA£) = ex P 
Z2AO = ex P 



a 



6 f + a*b 



exp 



b ] - ia*b 



£(l-T)|a|V_(s) 



exp 



?(T)|a|V + W 



where we defined two auxiliary functions 

a±(s) = cosh 2s ± cos(2# — (p) sinh 2s 



(23) 
(24) 

(25) 



and denoted by the phase of the squeezed vacuum. In the first factor of the generating 
operators we recognize the intrinsic quadrature operator 



*9 



6t + 



y/2 



(26) 



with the rotated quadrature xe+n/2 in Z 2 ^{0- The second factor of the generating operators 
may be seen as an operator ordering parameter introduced by Cahill and Glauber | i~9fl . By 
varying the filter parameters we can select a particular ordering, for example for 



0±(S 



1 + <t±(s) 

the creation and the annihilation operators are ordered antinormally 
Z lT (0 = exp {{i/2)a*VTa) exp ((e/2)av / Ta t ). 



(27) 



(28) 



For the discussed experimental NFM scheme, no selection of the parameters can lead to a 
normal ordering. 

The operational quadratures might be calculated explicitly 



y-(n) 



X. 



(n) 
2.F 



2z 
f 



(J- 18) 



H, 



T 



l-T ajs) 



2% 



a+(s) H n \ i 



1 — T Xo+tt/2 

~r~^(s) 



(29) 
(30) 



where H n (z) is the n-th Hermite polynomial and we have scaled the results by setting 
a = a/2. For the purpose of further discussion we write explicitly the first two operational 
moments for both quadratures 



y-(i) 



Y(2) 



Xq 



T 



x 



a_{s) 

2 l-T 



X. 



(i) 

2T 



(T-[S) 2 
x 9+tt/2 



«7_(s), 



<X+(s) 



X« = (1 - T) 



.1; 



+ 7T/2 

o-+(s) 
7 



T 



(31) 

(32) 
(33) 

(34) 



We see that for T = 1 or T = we are restricted to the measurement of either x[ n Jr or X^jr, 
which are, in this case, equal to the powers of the intrinsic operators. 

For simplicity, let assume that the phase of the squeezing vacuum is such that 29 + (ft = 0, 
then the functions cr±(s) becomes simpler 

°±(s)=e ±2s , (35) 



and the operational spread defined as 5Xijr = \J (X- 2 Jr) — (X- 1 ^) 2 reads for both quadratures 



e -2s 



{5X lT f = Te 2s (Ax e ) 2 + (1 — T) (36) 



{5X 2r f = (1 - T)e~ 2s (Ax e+7T/2 ) 2 + T ^, (37) 

where Axe is the intrinsic quantum dispersion. 

We note that for the selected filter parameters one of the operational quadratures is 
squeezed and the second one is enhanced. This behavior is the same as for the intrinsic 



quadratures |1| , however the measuring device introduces an additional factor contributing 
to the whole operational spread. Each operational spread contains two parts, the first part 
is the intrinsic uncertainty weighted by the factor T (or 1 — T) and the second part is the 
noise with the weight 1 — T (or T) introduced by the experimental device. The origin of 
these additional terms is clear. The NFM experimental setup enables us to perform a joint 
measurement of two quadratures in the extended Hilbert space. In the reduced space of 
the signal mode we pick up an additional noise due to the leaking noise mode through the 
unused ports in the NFM apparatus. 

By changing the transmittance coefficient we can select which quadrature will be mea- 
sured and how precise this measurement will be. If we set T = or T = 1 we measure only 
one quadrature and the moments of the operational operators equal to the powers of the 
intrinsic homodyne observables. The operational result reduces to the intrinsic observable 
in this case because we have assumed perfect photodetectors with 100% efficiency. 

The non ideal detectors can be described by a homodyne setup with perfect detectors, 
and a beam splitter with a certain effective transmittivity. The beam splitter mixes the 
incoming signal with the vacuum field and the results of the calculations are formally the 



same as those carried on for the non ideal detector |14|,|20]. 

It's worth to notice that in fact the NFM experimental scheme may be used for such a 
description of a quadrature measurement with a non ideal photodetector. In such a case we 
may forget about the upper right part of the NFM apparatus and see that the setup reduces 
to the situation described in In such a scheme there is a beam splitter in front of a pair 



of ideal detectors. If we put the squeezing parameter equal to 0, the operational observables 
obtained in this paper reduce to those obtained in 

We see that there are two possible and mathematically equivalent descriptions of the 
operational measurement performed with help of the NFM device. In one description the 
experiment is a simultaneous measurement of two operational quadratures with the help of 
ideal photodetectors with the NFM apparatus noise. In the other description, we perform a 
standard homodyne experiment, with a non unit quantum efficiency of the photo-detectors 
simulated by the beam splitter. Such a beam splitter is mixing the probe signal with the 



S 



vacuum state and the effective transmittance coefficient plays the role of the photodetector 
efficiency. 

Another possible choice of the parameters leads us to the operational momentum and 
position operators derived in ||, indeed for T = ~, s = and |a| = 2, 9 = our results can 
be reduced to the special case of the results from this reference, with the filter being in the 
vacuum state. 

The of operational operators x[ n Jr become then the operational moments of the position 
operator whereas X^p the operational moments of the momentum operator. 



V. SQUEEZED QUANTUM TRIGONOMETRY 

Originally the NFM experiments have been devoted to the measurement of the relative 
phase between a high-intensity classical field and a low-intensity laser field. The quantities 
which were investigated in these experiments were built out of the field quadratures fll9D 
and ( pUP ( with transmittance T = |). These operators have been constructed in such a way 
that they can play the role of the commuting "sine" and "cosine" operators. As a result the 
phase of the quantum probe field has been defined operationally, avoiding all the problems 
associated with the construction of a Hermitian phase operator in quantum mechanics [JJT] . 



It was shown |yj that one can associate with the NFM experiment an operational phase 
operator. 

In this section we extend this theoretical analysis for the case of the noise field being in 
a squeezed vacuum state (the original NFM experiments have been done with the vacuum 
state). We will show how to define an operator corresponding to any function of the phase. 
As a example we derive operators corresponding to the trigonometric functions. These 
operators form the squeezed quantum trigonometry of operational observables for the NMF 
setup. 

We shall start our discussion, introducing an amplitude marginal of the POVM for the 
NFM device 

fi((p) = J dlT{a = \fle i<p ). (38) 
From this POVM, we construct the corresponding phase propensity 

Prfa) = (39) 
which is a periodic function of phase, normalized to unity 

/^PrfoO = l. (40) 



A. Phasor Basis 

We define here the set of operational operators, which in the classical limit can be reduced 
to the Fourier expansion basis. The problem of constructing the operators corresponding to 
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the classical phase dependent quantities has been solved in Ref. Because every periodic 
function can be expressed in terms of its Fourier decomposition 



/(¥>)= E Ac**, ( 41 ) 

fc= — CO 

we will construct the operational quantum analogs of e lkv — so called phasors — in the 
following way 

7^={E {k) ), k = ±l,±2,.... (42) 
The reality of the propensity implies 

£(-fc) = £(m (43) 

because of the normalization condition (f40 ) we have = 1. 

As we have seen in the previous sections, the function f(ip) corresponds in the NFM 
device to the following operational operator 

CO 

F= Y: hE {k \ (44) 

k=— co 

which is a function of the boson creation and annihilation operators of the probe field. 

Because we have derived the POVM operator T for the NFM setup with an arbitrary 
noise field (fTIj), we can derive the explicit expression for the phasor E^ straight from 
the definition ([J). Once these operators are known, the operational quantum operator 
can be calculated using a Fourier series expansion. Because of this we will call the set 
{E^ k \ k = 0,1,2,...} the phasor basis. The definition of the phasors introduced in f22 
requires, that these operators become in the classical limit functions of the phase only. 

In order to find the phasor basis we start with the formula obtained in JTTj] 

= (-) Tr v {u k S(s,<j>)\0 v )(0 v \St(s,<j>)}, (45) 

where we have replaced the vacuum state of the filter by the squeezed vacuum state. 

The complex number a denotes the amplitude of the classical reference field of the local 
oscillator. The unitary operator U is defined as follows 



„ X l+ iX 2 ^X{ + X{ 



where Xi are given by ( |i~9D and (^) with the beam splitter transmittivity T — ~. Simple 
algebra gives 



V 1 / — — + V 1 



U = vW« r , 5 = V a */« n ■ (47) 
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We recognize in this expression the unimodular operator (H). In this way we obtain an 
equivalent expression for & k > 



Tr, 



(b + v 



(b + vt) (6+ + v) 



r S(s,<P)\0)(0\S\s,<P) 



(48) 



with k = 1, 2, 3, Negative values of in (£||) may be derived using fl43"|) . 

Because operator U obeys the assumptions imposed on the operator A(b + v^, ¥ + v) 
defined in (^) we can write the formula for the phasor using ([H]) as follows: 



d 2 a a k 



\a, s) (a, s\ 



(49) 



obviously F(s, <fi) = \a,s)(a,s\ is the POVM associated with the described experimental 
scheme. It is clear, that the POVM, contrary to spectral measure, depends on experimental 
device. In fact, for each value of squeezing parameter s we obtain a different propensity a 
different POVM and a different phasor basis, even though the probe field remains unchanged. 
An exact formula for the phasor may be derived straight from (|49"1) , 



S(s, -(, 



{b cosh s — we % ^ sinh s) k 



fet cosh s — be~ ir t' sinh s) 



(50) 



where : denotes the antinormal ordering. When the squeezing parameter s tends to zero the 
above formula reduces to the one from [0 . 



B. Properties of the Squeezed Quantum Trigonometry 

In this section we find various phase space representations for the simple combinations 
of the phasor basis. Using these relations we derive the trigonometric operational operators. 

These trigonometric operators will be defined in accordance with fl44j). We have the 
following operators 

5(1) = I(£(D _ £(-i)) , £<2) = I _ I(£(2) + £HQ) ; 

(7(D = + £;(-!)) , c® = ± + ±(E&> + &-*'>), (51) 

that correspond to the trigonometric functions sin <p, cos (p, sin 2 p> and cos 2 ip respectively. As 
it was stressed earlier and follows from definition of the phasors (C^ 1 )) 2 ^ C^. Consequently, 
the Pythagorean theorem sin 2 <p + cos 2 p = 1 is obeyed by S^ 2 \ C^, but not by (S^) 2 and 
(C^) 2 ). 

In order to investigate the properties of the phasors we evaluate their phase space repre- 
sentations corresponding to a certain arbitrary 77-ordering of boson operators, such a phase 
representation may be defined as follows [19] 
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m*Tr {& k \s,<l>)b{p)}. 



(52) 



The number r\ lies between minus infinity and plus one and P^ denotes the correspond- 
ing quasidistribution. Especially interesting cases are p( _1 ) ; p(°) ; p( -1 ) ; which correspond 
to the Q-representation (normal ordering), the Wigner function (symmetric ordering) and 
the P- representation (antinormal ordering), respectively. In this formula the phase space 
integration can be done and result we obtain 

d k X 



n— h. 



oo V7T 



x exp \ —X 2 Q 



E|/?| 2 + -A 2 sinh 2s (e^/3 2 + e^P* 2 ) 



where 



n 



l-A 



2^+1 



and 



A cosh s 



A I cosh s — 



l-A 



2^ + 1 



+ A cosh s 



77 + I 



+ 1. 



(53) 



(54) 



(55) 



This result is meaningful only if Q is real, otherwise one of the integrals, that appeared 
earlier, is divergent. Because the most interesting cases are the Q, P-representations and 
the Wigner function we have 

-l 



n 



(A 4 + 2A 2 ) cosh 2 s + 1 



iA 4 + A 2 (2sinh 2 s + 1) + 1 
1 - (A 4 + 2A 2 ) sinh 2 s 



for rj = — 1 
for rj 
for 7] 




1. 



(56) 



For r\ — 1, n becomes a complex number, so the P-representation for the squeezed phasors 
does not exists unless we tend with s to zero. In this limit we retrieve the original NFM 
phasors from O]. In our further analysis we shall restrict the calculations to the Wigner 
function only, but all the results may be easily rewritten for the respective Q-representation. 

If we write (3 = y/le 16 and use (|53|), the definitions fl5"T|) give for rj = the "sine" and the 
"cosine" Wigner functions 

-A 2 /[l - |A 2 + A 2 cosh 2 8 (1 + tanh s cos(0 + 26))} 1 



dX 



exp 



iA 4 + A 2 (l + 2sinh" s) + 1 



x- 



\fl cos 9 — ^A 2 cos 9 + A 2 cosh 2 s(cos 9 + tanh s cos(0 + 9)) 



and 



00 dX 



7T 



exp 



[iA 4 + A 2 (l + 2sinh 2 s) + l]f 

-A 2 /[l - f A 2 + A 2 cosh 2 s(l + tanh s cos(0 + 29))} 1 
iA 4 + A 2 (l + 2sinh 2 s) + 1 J 



(57) 



v^T sin ^ — | A 2 sin 9 + A 2 cosh 2 s(sin — tanh s sin(0 + 9)) 



[iA 4 + A 2 (l + 2sinh 2 s) + l] 



(58) 
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Similar expressions for the Wigner functions Syy (s, 0) and C^/(s, 0) can be derived from 
the phasor basis. Due to the lengthy mathematical form, we shall present in this paper 
only plots of these functions. The Wigner functions of the operational quantum trigonom- 
etry provide a simple framework for the physical discussion of the obtained results. In 
Fig. 2, we have plotted the Wigner function S^)(2, |). For field intensities tending to in- 
finity we should reproduce the classical trigonometry. In order to see this we recall that 
lim/^oo C\f^j e' 1 ? 2 ^ = 5(f(z)). After some simple algebra we see that /(A) in (|57|) and 
(g8|) has only one zero point at Ao = and the modulus of a derivative of /(A) in that point 
is one. As a result, using the well known relation between S(f(z)) and S(z), we see that 
in order to perform the classical limit, we have to put only A = in the integrands of the 
Wigner functions. In this limit 

7(1) 



r lim S { Q > olW (s,(f)) = sin6 
lim Cq ) ovW (s, 

I— >oo 



lim C£l TW (s,(j)) = cos6. (59) 



This result is also seen on Fig 2. where indeed for large values / we reconstruct the classical 
trigonometric function. 

As it might have been expected, in the classical limit the dependence on the squeezing 
parameters disappears, in agreement with the fact that squeezing is a purely quantum effect. 

Similarly we can prove, that S^(s,(P) and CW(s,<j>) posses a proper classical limit. For 
these functions we present only graphs. As it can be seen on Fig. 3 and Fig. 4, the classical 
limit of the Wigner function of C^ 2 '(s, 0) reproduces the cos 2 ip. 

We shall discuss now the limit of a very small probe field intensity. Let assume for a 
while, that the phase of squeezing parameter is zero, and I — ► 0. In this case we can 
expand the exponential factor of the integrand in (|58| ) and (|5TD. Keeping only the first term 
we obtain 

C$(s,0) 7 =° y/lA+(s,(/) = 0)cosfl 

S${s,(f>) *= VlA-{s,<f> = 0)sin6 (60) 

where 

* , i ^ f°° dX 1 + (l/2)A 2 e ±2s 

A±(s,j> = 0) = -=- K - J —^ = 3 61 

J-oo v^F [lA 4 + A 2 (l + 2sinh 2 s) + 1]§ 

are two different amplitudes for the "sine" and the "cosine" Wigner functions. These ampli- 
tudes depend only on the squeezing parameter. This result shows the effect of the amplitude 
squeezing near the vacuum. From (|61|) we see that for = S^(s, 0) is strongly squeezed 



then Cyy(s, 0) (i.e. the amplitude of sine is smaller then that for cosine). Obviously for 
1 = the amplitudes are exactly equal to zero, which is in agreement with the common 
intuition that the phase of a light field in the vacuum state is randomly distributed. 

Changing the squeezing phase to = 7r, we reverse the amplitude squeezing from the 
sine to the cosine function. It is simple to verify, that the amplitudes for different values 
of are connected in the following way A + (s, = 0) = A-(s, = tt) and A~(s, = 0) = 
A + (s, = 71"), we see that in this case, that Cy){s, 0) is strongly squeezed then S^(s, 0). 
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For arbitrary values of and s the separation of the cosine Wigner function into an 
amplitude, and a purely angle dependent part is no longer possible, but a clear squeezing of 
the amplitude can also observed comparing, for example, Fig. 3 (s=0.5) and Fig. 4 (s=1.5). 
For s = 2 (Fig. 2) the operational phasor is by far more squeezed then the two functions 
form Fig. 3 and Fig. 4. Note that the range of the intensity / is almost twice bigger then 
the range for Fig. 3 and 4. 

Similarly we can find an asymptotic expression for Cyy (s, 0). In the limit of small / 

C^ ) (s, ( j>) = ^(l-c(s, ( j>)), (62) 

where c(s, 0) is an /-independent function of the squeezing parameter 

, „ rfAirfA 2 (l/2)(l + (l/4)A 4 + A 2 (l + 2sinh 2 s))sinh2scos0 , , 

c(s,0) = / § , (63) 

J-oo n [iA 4 + A 2 (l + 2sinh 2 s) + l]5 

where A is understood here as a norm of the two component vector A = [Ai, A2]. It's easy 
to check, that c(s, 0) tends either to unity (0 = 0) or to minus unity (0 = ir). As a result, 
for small intensities C^(s, 0) becomes zero (0 = 0) or one (0 = it). For = tt/2, 3tt/2 we 
have c(s, 0) = and lim/_,oC^(s,^) = 1/2. 

For small / the squeezing influences the system very strongly. If the squeezed phase 
equals to 7r/2, the Wigner function tends to 1/2 (Fig. 3), whereas for = it takes 
values near zero (Fig. 4). Such a dramatic change of the cosine quadratures occurs because 
in the limit of small J, purely quantum effects of the squeezed vacuum are important. The 
squeezing allows one of the quadratures to be reduced below the vacuum level represented by 
a uniformly distributed random phase. The uniform distribution of the phase corresponding 
to the vacuum state leads to an operational Wigner function for / = equal to ~. For a 
squeezed vacuum, this uniform random-phase distribution is modified |23j and a significant 
change of the operational quadrature is possible. In fact fluctuations below | in the 
Wigner function exhibit the quantum nature of the squeezed vacuum. Those results are even 
more clear in view of our previous analysis of operational quadratures (|l^) and (j20l) . We see 
straight from the definition of U ((|§), that operational operators C™(s, 0) and S^ 2 \s, 0) are 
just normalized homodyne quadratures, so their dependence on squeezing phase and s may 
be understood as if they were homodyne quadratures. These operators are also quantum 
analogs of certain classical functions of the phase of the electromagnetic field. Because of 
this they posses a much richer structure and this enable, for example, to investigate the 
phase properties of the quantum optical fields. 



VI. CONCLUSIONS 



We have presented an operational theory of the eight-port homodyne detection scheme. 
For such a measuring apparatus, we have derived the operational quadrature operators for 
an arbitrary noise field leaking through the unused port. In the framework of the opera- 
tional approach we have derived the quantum propensity, the POVM and the corresponding 
operational observables for the NFM device. For the noise field being the squeezed vacuum 
we have constructed the quantum trigonometry and the corresponding Wigner functions. 
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FIGURES 




FIG. 1. The NFM apparatus. The photodetectors are labeled by numbers 3, 4 and 5, 6, whereas 
the beam splitters are represented by dark thickened lines 




FIG. 2. Plot of the Wigner function of the operational operator S^(s = 2, <p = tt/2). 
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FIG. 3. Plot of the Wigner function of the operational operator C^(s = 0.5, <p = tt/2). 




